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Abstract. Within the next few years, GAIA and several instruments aiming at imag- 
ing extrasolar planets will see first light. In parallel, low mass planets are being searched 
around red dwarfs which offer more favourable conditions, both for radial velocity de- 
tection and transit studies, than solar-type stars. Authors of the model atmosphere code 
which has allowed the detection of water vapour in the atmosphere of Hot Jupiters re- 
view recent advancement in modelling the stellar to substellar transition. The revised 
solar oxygen abundances and cloud model allow for the first time to reproduce the pho- 
tometric and spectroscopic properties of this transition. Also presented are highlight 
results of a model atmosphere grid for stars, brown dwarfs and extrasolar planets. 



1. Introduction 

Since spectroscopic observations of very low mass stars (late 80s), brown dwarfs (mid 
90s) and extrasolar planets (mid 2000s) are available, one of the most important chal- 
lenges in modelling their atmospheres and spectroscopic properties lies in high temper- 
ature molecular opacities and cloud formation. K dwarfs show the onset of formation 
metal hydrides (starting around T sS ~ 4500 K), TiO and CO (below T eS ~ 4000 K), 
while water vapour forms in early M dwarfs (T e s ~ 3900 - 2000 K), and methane, am- 
monia and carbon dioxide are detected in late-type brown dwarfs (T e g ~ 300 - 1600 K) 
and in extrasolar giant planets. The latter are either observed by transit (r e ff ~ 1000 - 
2000 K depending on the spectral type of the central star and the distance to the star) or 
by imaging (young planets of T e g ~ 300 - 2000 K depending on their mass and age). 

The modelling of the atmospheres of very low mass stars (hereafter VLMs) has 
evolved (as here illustrated with the development of the PHOENIX atmosphere code, 
which has allowed the detection of water vapour in extrasolar planets' atmospheres by 
Barman et al. 2007, 2008) with the extension of computing capacities from a n ana- 
lytica l treatment of the transfer equation using moments of the rad iation field (Aliardl 
1990), to a line-by-line opacity sampling in spherical symmetry (I Allard etaL 1 199"/ 



Hauschil dtetaDl 1 999alfbh . and more recently to 3D radiation transfer (ISeelmann et all 



2010). In parallel to detailed radiative transfer in an assumed static environment, hy- 
drodynamical simulations have been developed to reach a realistic representation of the 
granu lation and its induced line shifts for the sun and sun-like stars (see iFreytag et al.l 
l201lh by using a non-grey (multi-group binning of opacities) radiative transfer with a 
pure blackbody source function (scattering is neglected). 
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Figure 1. Synthetic spectra compared to the IR SED of VB10 using iden- 
tical model parameters (T e ff = 2800K, logg = 5.0) and a resolution of 50 A 
with d ifferent water vapour o pacity sources: the Base grid by Alla rd & Hausc hildt 
(119951) u sing iLudwid (11971b: a test using the 1994 version of th e Niel s Bohr 
Insti tute dJ0rgensen et all 1200 ll) : the NextGen grid by lAllard et all dl994 Il997h 
and [H auschildt et al. ( 1999a b) using the University Colle ge London d atabase 
(Schr vberetalJ 11995b : and the AMES-Cond/Dusty grid b v lAllard et all (12000b: 
Allar d et al.l (12001) using the NASA-Ames Center database (Partridge & Schwenke 
19971) . All models (except the NextGen/UCL case) underestimate the flux at K (ca. 



2.0 - 2.4 ftm) by 0.1 to 0.2 dex. 



2. Molecular opacities 



While earlier work has been developed for the study of red giant sta rs, the pioneer - 
ing work on t he m odel l ing of VLM atmospheres has been provided by iMould d 19751) . 
Allard d!990h and lKul (Il99ll) using a b and model or Just Overlapping Line Approxi- 
mation (J OLA) opacities developed by iKivel et all d 19521) and adapted for astrophysi- 
cal use bv lGoldenldl967l ). More realistic model atmospheres and synthetic spectra for 
VLMs, brown dwarfs and extrasolar planets have been made possible thanks to the de- 
velopment of accurate opacities calculated often ab initio for atmospheric layers where 
temperatures can reach 3000 K. The process of improvements was especially remark- 
able in the case of water vapour line lists. Indeed, water vapour has seen an important 
evolution through the years from band model approximations to straight means based 
on hot flames experiments, and then to ab initio computations. Nevertheless, the at- 
mosphere models have failed to reproduce the strength of the water bands that shape 
the low resolution (R < 300) infrared spectral energy distributions of M dwarfs. At the 
lower temperatures of brown dwarfs methane and ammonia rival the effect of water. 
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The discrepancies in the model synthetic spectra were therefore believed to be due to 
inaccurate or incomplete molecular opacities. In particular water vapour was suspected 
because the discrepancies were observed at infrared wavelengths in the relative bright- 
nesses of the flux peaks between water vapour bands. As can be seen from Fig.[T]where 
the models are compared to the infrared spectrum of the M8e dwarf VB10, the water 
vapour opacity profile which shape this part of the spectrum has strongly changed over 
time with the improvement of computational capacities and a better knowledge of the 
interaction potential surface. And the most recent ab initio results confirm the earliest 
hot flames laboratory experiment results by lLudwigl ( 197 lh . But in general, most opac- 
ity profiles produce an excess opacity (or lack of flux in the model) in the K bandpass. 
Only the UCL1994 line list (due to incompleteness, and with much of its deviations 
cancelling out over the bandpasses) could produce seemingly correct J - K s colours. 



3. The revised solar abundances 

Model atmospheres for VLMs and in general for other stars assume scaled solar abun- 
dances for all heavy elements, with some enrichment of a-process elements (the result 
of a "pollution" of the star-forming gas by the explosion of a supernova) when appro- 
priate in the case of metal-depleted subdwarfs of the Galactic thick disk, halo and glob- 
ular clusters. The revision of the solar abundances based on radiation hydrodynamical 
simulations of the solar atmosphere, on improvements in the quality of the spectro- 
scopic observations of the Sun, and in its detailed line profile analysis by two sepa- 



rate groups using in dependent hydro codes and spectral synthesis codes dAsplund et al. 



120091 : ICaffauet all 1201 lh yield an oxygen red uction of 0.11-0.19 d ex (up to 34%). 
compared to the previously used abundances of iGrevesse et al. (1993). Since the over- 



all SED of late K dwarfs, M dwarfs, brown dwarfs, and exoplanets is governed by 
oxygen compounds (TiO, VO in the optical and water vapour and CO in the infrared), 
the input elemental oxygen abundance used in the equation of state is of major impor- 
tance. Fig.|2] shows an example of these effects for the optical and infrared SED of 
the M5.5 dwarf system Gl 866. However at other effective temperatures even stronger 
photometric effects can be seen, where the near-IR SED of different models diverges 
more (see Fig.[3]>. The comparison shows significant improvement compared to older 
models shown in Fig.Q] except for excess flux in the H bandpass near 1.7 fim due to 
incomplete FeH opacity data for this region. The comparison has particularly improved 
in the Wing Ford band of FeH near 0.99 /urn, and in the VO bands thanks to line lists 
provided by B. Plez (GRAAL, Montpelier, France), although inaccurate or incomplete 
opacities are st ill affecting the models at optical wavelengths (e.g. the TiO line list by 
lLanghoff|[l997h . 

Fig. [3] compares th e theoretical isochrones (assuming an age of 5 Gyrs) to the 
Casagrande et al. (2008) r e .ff e stimates and reveals that the NextGen models (lAllard et all 



1997NHauschildt et all 19993 15) systematically and incre asingly overestima te T e s through 



the lower main sequence, while the AMES-Cond/Dusty dAllard et al.ll200lh models on 



the contrary underestimate T e ff as a function of / - K s colour. This situation is relieved 
when using the current models (labeled BT-Settl in the figure) based on the revised 
solar abundances, and the models now agree fairly well with most of the empirical es- 
timations of r e ff. The current model atmospheres have not yet been used as surface 
boundary condition to interior and evolution calculations, and simply pr ovide the syn- 



thetic colour tables interpolated on the published theoretical isochrones dBaraffe et al. 
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Figure 2. A BT-S ettl synthetic spectrum with log g-5.0, and solar metallicity by 
Aspl und et al.l (120091) is shown as lig ht sohd line, compared to the combined SEE) 
of the red dwarf triple system Gl 866 dLeinert et al.lll990tlDelfosse et al.ll 1999b . The 
observations were combined from a Mt. Stromlo optical spectru m (M. Bessell, priv . 
comm.) and SpeX infrared spectrum taken at the NASA IRTF (Ra yner et al.l l2009) 
(thick cu rve). For comparison a model using the same parameters and physical setup 
with the iGrevesse et alj (1 19931) abundances is shown as a dashed line. The mod- 
els have been scaled to the observed absolute flux assuming two equal r e ff=2920K 
components of 0. 157 solar radii and a third with T e ff =2700 K and 0. 126 R . 



998). Even if the atmosp heres partly control the cooling and evolution of M dwarfs 



(IChabrier & Baraffd ll997h . differences introduced in the surface boundary conditions 



by changes in the model atmosphere composition have negligible effect. 



4. Cloud formation 



One of the most important challenges i n mod elling these atmospheres (below 2600K) 
is the formation of clouds. iTsuji et al.1 (I1996T) had identified dust formation by recog- 
nising the condensation temperatures of hot dust grains (enstatite, forsterite, corun- 
dum: MgSi03, Mg2Si04, and AI2O3 crystals) to occur in the line-forming layers 
(t w 10~ 4 - 10~ 2 ) of their atmospheres. The cloud composition, according to equi- 
librium chemistry, is going from zirconium oxide (Zr02), to refractory ceramics (per- 
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Figure 3. 




for M dwarfs by [ Casagran de et all (|2P08) and brown 
dwarfs bv lGolimowski e^ alJ(l2004l):IVrba et al. | (120041) are compared to the NextGen 
isochrones for 5 Gyrs (iBaraffe et al.lll997Lll998h using various generations of model 
atmospheres : Nex tGen (thick line), the limiting case AMES-Cond/Dusty cases by 
Allard et alj d2001[) ( dotted and dashed lines), the current BT-Settl models using the 
Asplund et al.l ([2009) solar abundances (full line). The Gl 866 system fitted in Fig. [2] 
is highlighted by darker colours and shown with its relatively large photometric error 
bars at / - K« = 0.9. 



ovskite and corundum; CaTi03, AI2O3), to silicates (e.g. forsterite; Mg2Si04), to 
salts (CsCl, RbCl, NaCl), and finally to ices (H 2 0, NH 3 , NH4SH ) as brown dwarfs 



cool down over time fro m M through L, T, and Y spectral types (|Allard et al. 2001 



lLodders & Fegleyll2006l ). This assumed (by Allard et al. 2001) sub-micron-sized crys- 
tal formation causes the weakening and vanishing of TiO and VO molecular bands (via 
CaTi03, Ti02, and VO2 grains) from the optical spectra of late M and L dwarfs, reveal- 
ing CrH and FeH bands otherwise hidden by the molecular pseudo-continuum, and the 
resonance doublets of alkali transitions which are only condensing onto salts in late-T 
dwarfs. The scattering effects of this fine dust is Rayleigh scattering which provides 
veiling to the optical SED of late-M and L dwarfs, while the greenhouse effect due to 
the dust cloud causes their infrared colours to become extremely red compared to those 
of hotter low mass stars. The upper atmosphere, above the cloud layers, is depleted 
from condensible material and significantly cooled down by the reduced or missing 
pseudo-continuum opacities. 

One common approach has been to explore the limiting properties of cloud forma- 
tion. One limit is the case where sedi mentation o r gravitational settling is assumed to be 
fully effici ent such as the case B of lTsujil (2002), the AMES-Cond o r condensed phase 
models of lAllard et ail (1200 ll). the clea r case of lAckerman & Marleyl d200ll) and the 
cloud-free case of Burrows et al. (2006). The other limit is the case where gravitational 
settling is assumed inefficient and d ust, often on ly forsterite, forms in equilibrium with 
the gas phase such as the case A of iTsuiil (120021). the AMES-Dusty or dusty models of 



Allard et all (|2001 ), the cloudy case of lAckerman & Marleyl (l200ll) . or the case B of 
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Burr ows et al. (2006). These limiting cases of maximum dust content agree in describ- 



ing the evolution of brown dwarfs from a molecular opacity governed SED towards a 
blackbody SED below 1500K. This description was suitable, at least in the case of the 
AMES -Dusty models, in reproducing the infrared colours of L dwarfs. The cloud-free 
limiting case on the other hand allowed to reproduce to some degree the colours of 
T dwarfs. Fig.B] show s this situation for the AMES-Cond/Dusty limiting case models 
of lAllard et all ( 20011) compared with the effective temperatures estimates obtained by 
integration of the observed SED dGolimowski et al.ll2004l : IVrba et alJl2004h . 

The purpose of a cloud model is therefore to go beyond these limiting cases and 
define the number density and size distribution of condensates as a function of depth 
in the atmosphere. The discovery of dust clouds in M dwarfs and brown dwarfs has 
therefore triggered the development of cloud model s building up o n pioneering wo rk in 
the context of plane tary atmospheres developed by iLewisI (1 19691) . IRossowI ( 1978 ). and 
Limine et al. (1989). The Lewis model is an updraft model (considering that condensa- 



tion occurs in a gas bubble that is advected from deeper layers). By lack of knowledge 
of the velocity field and diffusion coefficient of condensates in the atmospheres of the 
planets of the solar system, Lewis simply assumed that the advection velocity is equal 
to the sedimentation velocity, thereby preserving condensible material in the conden- 
sation layers. This cloud model did not account for grain sizes. Rossow on the other 
hand developed characteristic timescales as a function of particle size for the main 
microphysical processes involved (condensation, coagulation, coalescence, and sedi- 
mentation). The intersections of these characteristic timescales gives an estimate of the 
condensate number densities and mean grain sizes. However, this model made several 
expli cit assumptions conc erning the efficiency of supersaturation, the coagulation, etc. 

iHelling et all d2008al) have compared different cloud models and their impact on 
model atmospheres. Most cloud models define the cloud base as the eva poration laye r 
provided by the e quilibrium chemistry. In the unified cloud models of iTsujil ([2002); 
Tsuji et all (12004) a parametrization of the radial location of the cloud top by way of 



an adjustable parameter T cr jt was used. This choice permits to determine the cloud 
extension effects on the spectra of these objects but does not allow to reproduce the 
stellar-substellar transition with a unique value of T cr jt since the cloud extension de- 
pen ds on the atm o spher ic parameters. 

lAllard et all d2003l) using PHOENIX and the index of refraction of up to 40 con- 
densible species, have applied the Rossow cloud model, ignoring coalescence and co- 
agulation, and comparing the timescales of condensation, sedimentation and mixing 
(extrapolated from the convective velocities into the convectively stable layers), and 
assuming efficient nucleation (monomers equilibrium densities). The cloud model was 
then solved layer by layer inside out to account for the sequence of grain species for- 
mation as a function of cooling of the gas. But this version of the BT-Settl (with gravi- 
tational settling) models did not allow for the formation of enough dust in brown dwarf 
atm ospheres due to a too conserv ative prescribed supersaturation value. 

Ackerman & Marleyl (1200 ll ) have solved the particle diffusion problem of con- 



densates assuming a parametrized sedimentation efficiency / se d (constant through the 
atmosphere) and a mixing ass umed constant and fixed to its maximum value (maximum 
of the inner convection zone). ISaumon & M arlev (2008) found that their models could 
not produce the co l or chan ge with a single value of / se d. 

IHelling et all (l2008bl) use the PHOENIX code to compute the Drift-Phoenix mod- 
els. The cloud model used, in the contrary to all other cases mentioned, studies the 
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nucleation and growth of grains as they sediment down into the atmosphere. This cloud 
model determines the number density and size distribution of grains by ID nucleation 
simulations, and the resulting distribution is read in by PHOENIX which computes the 
resulting opacities and radiative transfer. These models solve the nucleation problem 
but only for the assumed monomer types and have been successfully applied to fit the 
dusty atmospheres of L dwarfs, but the reversal i n IR colours observed for the L/T 
transition could not be explained dWitte et al.ll2011 ). 

None of these models however treated the mixing properties of the atmosphere and 
the resulting diffusion mechanism realistically enoug h to reproduce the b rown dwarf 
spectral transition without changing cloud parameters. iFreytag et al.l ( 2010h have there- 
fore addressed the complementary though important issue of mixing and diffusion in 
these atmospheres by 2D radiation hydrodynamic (hereafter RHD) simulations, using 
the PHOENIX gas opacities in a multi-group opacity scheme, and forsterite with geo- 
metric cross-sections. These simulations assume efficient nucleation, using monomer 
densities estimated from the total available density of silicon (least abundant element 
in the solar composition involved in forsterite). They found that gravity waves play a 
decisive role in clouds formation, while around T e s < 2200 K the cloud layers become 
optically thick enough to initiate cloud convection, which participate in the mixing. 
Overshoot can also be important in the deepest layers. 

These RHD simulations allow an estimation of the diffusion processes bringing 
fresh condensible material from the hotter lower layers to the cloud forming layers. We 
have therefore updated our cloud model (BT-Settl models) to account for the mixing 
prescribed by the RHD simulations. Another important improvement concerns the su- 
persaturation which as been computed rather than using the fixed conservative value 
recommended by Rossow. One can see from Fig.@]that the late-type M and early-type 
L dwarfs behave as if dust is formed nearly in equilibrium with the gas phase with ex- 
tremely red colours in some agreement with the BT-Dusty models. The BT-Settl models 
reproduce the main sequence down to the L-type brown dwarf regime, subjected in the 
K bandpass to the greenhouse effect of dust clouds, before turning to the blue in the 
late-L and T dwarf regime as a result of methane formation in the K bandpass. This 
constitutes a major improvement over previous models, and is promising that we can 
reach in the near futur a full explanation of the stellar substellar transition. 

Diffusion has also been held responsible for deviations in ultracool atmospheres 
from gas phase chemical equilibrium (CE), as noted i n early observat i ons of T dwarfs 
show ing an excess of carbon monoxide absorption ( Noll et all 1 19971 : [Griffith & Yelld 
1999). More recently, carbon dioxide (Tsuj i et ai1l2011 ). which was not expected at 



such low tempe ratures, has been de tected. Similarly, ammonia has been shown to be 
underabundant dSaumon et al. 120061) . This is understood as the result of slowing down 
of crucial chemical reaction steps, so that some important molecules (CH4, NH3) would 
not have the time to form in equilibrium while undergoing mixing, w hereas others 
(CO, C O?, N2) remain at enhanced abundances. The RHD simulations of IFreytag et al. 



( 20101) have allowed to understand the underlying mixing processes, obviating the need 
to describe them with an additional free parameter. 



5. Applications to Exoplanet Science 



Several infrared integral field spectrographs combined with coronograph and adaptive 
optic instruments being developed are coming online before 2013 (SPHERE at the VLT, 
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Figure 4. Same plot as Fig.[3]but zooming out and extending into the brown dwarf 
region of the diagram. This region below 2500K is dominated by dust formation (es- 
sentially forsterite and other silicates). The limiting cases AMES-Cond/Dusty model 
atmospheres provide a description of the span in colours of the brown dwarfs in this 
diagram for a given age (here 5 Gyrs). The BT-Settl models succeed in explaining 
even the most extreme colours of brown dwarfs. 



the Gemini Planet Imager at Gemini south, Projectl640 at Mount Palomar, etc.). The E- 
ELT 41m telescope in Spain due around 2020 will also be very ideally suited for planet 
imaging. The models developed for VLMs and brown dwarfs are a unique opportunity, 
if they can explain the stellar-substellar transition, to provide a great support for the 
characterisation of imaged exoplanets. We have therefore developed the BT-Settl model 
atmosphere grid to encompass the parameter regime of these objects (surface gravity 
around log g=4.0, T eS below 2000K). 

These planets are typically found at several dozens of AU from the star, and since 
the observatio ns are done in the infr ared the non-irradiated models can even be used di- 
rectly. Indeed, iBarman et al.1 (1200 ll) have shown that the effects of impinging radiation 
from a star on the planetary atmosphere are Rayleigh scattering of the stellar light by 
H2 molecules (or clouds if present) at optical wavelengths (below 1 /mi for solar type 
stars), while the impact on the interior and evolution properties becomes negligible for 
orbital distances exceeding 0.1 AU. Nevertheless, we are developing for 2012 irradiated 
models and the capacity to compute them via the PHOENIX simulator (see below). 



6. Summary and Futur Prospects 



We report progress of the development of a new model atmosphere grid for stars, brown 
dwarfs and young plan ets, named BT-Sett l. It has been c omputed using the PH OENIX 
code updated for: i) thelBarber et al.l(l2006l) BT2 water, thefHomeier et all d2003h STDS 
methane, the lSharp & Burrows! (l2007h ammonia and the lTashkun et all (|2004F CDSD- 
1000 CO2 opacity line lists, ii) the solar abundances revised by Asplund et al.l (2009), 
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and iii) a cloud model accounting for more detailed supersaturation and RHD mixing. 
The grid is covering the whole range of stars to young planets 400 K < T e g < 70,000 K; 
-0.5 < logg < 5.5; and -4.0 < [M/H] < +0.5, including values of the a-element 
enhancement (supernovae enrichment of the star forming material) between +0.0 and 
+0.6. Models are available at the PHOENIX simulator website http : //phoenix . ens-lyon . f r/simulator/ 
and are in preparation for publication to serve among others the GAIA, MUSE and 
SPHERE/GPI/P1640 instruments to come online in the near futur. Corresponding evo- 
lution models are expected for 2012. 

We found the previously used NextGen models to systematically overestimate T e s 
below 3500 K by as much as 500 K. The water vapour opacity profile has converged 
with the most recent line lists reproducing laboratory results, but could not explain this 
discrepancy. The solution came instead from the revision of the solar abundances which 
changes the strengh of the water vapour absorption bands, and therefore allows the 
reproduction of the spectroscopic and photometric properties of M dwarfs as late as M6. 
Later-type M dwarfs are affected by dust formation and cloud modelling is important to 
understand their properties. We find that the Rossow cloud model allows, with revisions 
to the supersaturation and mixing, the reproduction of the stellar-substellar transition. 
A small offset persists however in the M-L transition. It is possible that all the current 
cloud models are not efficient enough in producing dust at the onset of cloud formation 
regime. Detailed nucleation studies could allow in the futur to resolve this issue. Other 
uncertainties affect the current cloud modelling such as the assumption of spherical 
non-porous grains while grains form as fractals in the laboratory. Constraining the 
models remains therefore very important. 

Beyond cloud modelling and molecular opacities, model atmospheres for these 
objects require reaction rates for the most abundant molecules and/or most important 
absorbers. Furthermore, these atmospheres are composed of molecular hydrogen which 
constitute the main source of collisions. Also needed are therefore collision rates (by 
H2) and corresponding damping constants for the broadening molecular lines. 

In order to say something about the spectral variability of VLMs, brown dwarfs 
and planets, 3D global or "star in a box" RHD simulations with rotation will be re- 
quired. This is our current project supported by the French "Agence Nationale de la 
Recherche" for the period 2010-2015. 
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